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15. We shall now apply the equations of Article (12) to 
the Lunar Theory, in the following manner. We shall first 
suppose that the moon moves in the plane of the ecliptic, and 
afterwards prove the results obtained on this hypothesis to be 
true when the inclination is taken into account; we shall then 
determine the motion of the plane of the orbit, the mean 
motions of the perigee and node to the third order, and the 
parallactic inequality. 

By Article (12), we have 


_ dR 
and by Article (3), we have 
d’r du ,,., dudh 
de ae a at 
Therefore (1) becomes 
d*u uw r dR i1drdR 


wt*"2 Pe Pana © 


and since . = 3 (2) becomes 


d(h*) _ ‘ 
<a r oe 


16. To determine R, we have, by Article (14), putting 
cos (7, 7’) = p, 
neti R-m (= -5p), 
en 


1 





50 On the Lunar Theory. 


ee 273 r 
and — = —(1—2pz+2")* wherez=- 
anc ; = ( pz + 2°) whe - 
1 f 3.5 Ine — YR 
“oi \(2pe—2*)+ (ape - 2+ “TT aa 2°)'&e... 
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i} 


, P - 1 2 , oP — 3p 3 &e 
P 5 z coop 


ro r 
Hence, retaining only the lowest power of z or =, we have 
r 


(observing that p or cos (7, 7’) = cos (0 - 0’), since we suppose 
the moon to move in the plane of the ecliptic ) 








R-™ 4 ll 1+83 cos 2 (0 - 0) 
r r 4 
and therefore dh = i i Past (0- 0) eee | |S 
dy r° 2 
iR mr 3sin2(0-0 
a i En 
dd , 2 


We shall expand the second members of (3) and (4) as far 
as the third order of small quantities, retaining, however, only 
those terms of the third order in which the coefficient of 6 is 
nearly equal to unity or zero ; and this we shall do by substi- 
tuting for 7, 7’, and 6’ their values in terms of @ to the first 
order, namely 


r=a{1l-ecos(@-w)} #7 =a' {1-e' cos(0'-w)}... (7), 
0 = m9 + B + 2e' sin (m0 + B-w') ...... (8). 


17. To determine in in this manner, we have, by (4) and (6), 
F 





d (h’) ar 7 
—— = n2(60-0); 
dé re )s 
or, by (7), putting a 4 =m, 
: ua 
dh’) 


_ 3aum’ {1 — 4e cos (8 —w) + 8e' cos (6 — w')} sin 2(0- 0’). 

In the terms of this, which are multiplied by e and e’, we 
may put for 6’ its first approximate value, namely, mO + 3. 
In the remaining term, if we put for 6’ its value (8), expand 
to the first power of ¢’, and then transform products of sines 
and cosines into sums, in the usual way, it is evident that 
sin 2 (6 — 6’) will contain no term in which the coefficient of 
@ is nearly equal to unity or zero. It follows therefore that 
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we may suppose 6! to be equal to m@ + (3, simply, in all the 
terms of this equation. 

This being the case, the product cos (6 — w’) sin 2 (0 - 6’), 
transformed into a sum, gives no term to be retained, and the 
product cos (8 ~ w) sin 2 (- 6’) gives only } sin (0 - 20' + w): 
we have therefore 

A 3aum’ {sin 2 (0 - 6’) - 2e sin (0 - 20' + w)}. 
do 

Therefore, integrating, supplying the proper constant 
{au (1 - e*)}, and neglecting powers and products of m and e 
above the second order, we have 


3” 


h? = ay {1 -e@4+ > cos 2 (0 — 0) - 6m’e cos (0 - 26’ + ») - 


1 [ » 8m ' 
iy +e cos 2 (0 — 0')+ 6m’e cos (0 - 26 val. 
a | 2 J 


2 
18. To find - A eR , we have, by (5) and (7), 
h*® dr si 
m'r’ 1+ 3 cos 2(0- 0) 
apr” 2 





° ' » ny 148c0s2(0-0 
“aie {1-8ecos(9-w)+3e' cos(0'-w')}. en. 
a 
For exactly the same reasons as before, we may here suppose 
6 to be equal to m@ + 3 simply; and then, retaining only the 
proper terms, we have 
2 2 
“a cB = {1 - $008 2 (0-8) +$ ¢ 008 (0 —w) 
i j a Y ’ U U 
’ + $e cos (0 — 20 + w) - 3 e' cos (0 - w)}. 
1 dr dR ; : 
9. Lastly — — -——,we have, by (6) and (7), pro- 
19. Lastly, to find Fo do’ , by (6) )»P 
ceeding exactly as before, 
1 dr dR 3m'r dr 


it d0 dO apr” dO 


<< es (0 — w) sin 2 (0 - 6’) 
2a 


sin 2 (0 - 6’) 


2 
=— = e€ cos (0 - 27 + w). 
4a 


20. Let us now, for brevity, put 
2(0-@)=9, 0-27 +u=y, W-w = x) 
6’ here being simply m8 + B; 
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Also, for e cos (8 - w) put aw- 1; and then we have 
2 2 
ew lite mis. 
= =— +— i— — COs @ + 6e COS 
h’ a { 2 % v} 
r dR 3m m f 3 Ba. } 
—- = — =— Uu+ — 4-2-—Ccosgi—ecosy — — é COSyx> 
h* dé 2 a \ 2 ¢ 4 ¥ 2 xy 
1dr dR m 3 \ 
— —— ——=— 4—~—€COsy>. 
h' dd d@ a { 4 ¥ 
Hence, (3) becomes 
T’ 3m’ \ 
C wats _ 3m * 


dQ \ 2 / 


1 2 2 2 , . 
=- (1 +e — 2m’ - 3m’ cos d +  m’e cos ~ — 3 m’e' cos x.) 


Cc 


: 3m? iy, 
Assume now | putting 1 - . C 
\ 


ie {e cos (c0 - w) + A+ Boos ¢+C cos p+ D cos w} ; 
a 


and therefore, by (9), 


d’ - , 
— + Ue = — [c°.A + {c’ -(2 - 2my} Bos o 
a 


de? ; , ; 
+ {e?-(1-2m)*} Ccos py +(e’ - m’) D cos x]. 
Hence, equating coefficients, neglecting powers and pro- 
ducts of m and e above the second order, and putting 


1 
- (1+ & — 2m’) =u,, we have 
a 








2 15 sme 
uU=U, " +e cos (c8-w) + m’* cos p+ — cos p— 9 cosy} . 


21. To find ¢ in terms of 0, we have 


dt r 

dd ih 
in which we must substitute the values of A and r already 
found in Arts. (17) and (20), and then integrate. But since 
h is found by an integration with respect to 0, and ¢ by 
another integration with respect to 0, it is not sufficient to 
retain merely the terms of the third order in which the 


d(h’) 


ude bee 





coefficient of @ is small in the expression for , but we 


must also preserve all similar terms of the fourth order.* 





* In Airy’s Tracts and Pratt’s Mechanical Philosophy the terms of the third 
order only are retained, and therefore it appears to me that ¢ is found by an 
inconclusive method, so far as terms of the second order are concerned. 
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2 


v) 
-, we must 


dé 
substitute for r*, r* and 6’ their values to the second order, 
namely 





'To obtain these terms in the expression for 





, 5me 
real 140 4ecos(O-w)+ 5e*cos 2(0-w)— 4m’ cosg-- ony}, 
r=a' {1 — 3e' cos (6 — w)} ; 
(we may look upon e” as of the same order as m’). 
0 = m0 + B + 2e' sin (m0 + B - w') - 2me sin (0 - w). 


Now it is easy to see that, on account of the expression for 


a being multiplied by sin 2(@— 6’), the expansions of 


r and @' will give rise to no terms in which the coefficient of 
0 is small ; we may therefore suppose r’ = a’ and 0’ = m0 + B 
simply. But in the expansion of 7*, the term 5e’ cos 2(0 — w) 
combined with sin 2 (0 — 6’) gives a term to be retained, 
namely - 3 e’ sin 2(0'-w). We have therefore 





d(h’ ; 
a o 3aum {... ....- 36 sin 2(0'- w)}, 
Keanu {..... . — 2» me cos 2(8' - w)} 
fre Cave sins . + P me’ cos 2 (0' - w)}. 


Introducing this additional term in Art. (20), we find 
w=u,{...... P me cos 2(0' - w)}. 

22. We may now proceed to find ¢. 

Represent the values of w and h’ by wu, (1+ U) and 
au (1 — e’)(1+ H); then we have, by (10), 
dt ‘ 1 
—=- uu h? o4 = SS 
do (*) u,V{an(1—-e 
The values of U and H (neglecting in H the term of the 
third order in which the coefficient of @ is not small) are 


ay 2U+3U?-3H)...(11). 


U = e cos (cO - w) + m’ cos 


15 5me 3m’e’ 15me’ 


os ~ - cos x + —7— cos 2(0'~ w), 





38m 15me’ 


oe 4 cos 2 (9 - 6); 
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U* evidently contains one term of the second order and two 
of the third, namely, 
é cos’ (cO -—w), 2m’e cos (cM - w) cos ¢, 
15me* 
and - - os (cO - w) cos y. 


The second of these gives no term to be retained, the third 

gives (putting ¢ = 1 and observing that ~ = 0 - 20+ w), 

15me* — 
“ cos 2 (6 - w), and the first is equal to 
Le {1.4 cos 2(cO - w)}. 

Hence the coefficient of ~ term 15m*e cos 2(0' -— @) in 

(11) is—§+%+%, which = 0; this term therefore disappears. 

Therefore, collecting the abl of (11), and putting all the 





1 

constant part = —, we have 
n 

ndt 


dé 
as > ae 15me 
+ —- cos 2 (cB - w) - - cos p — 


= 1 - 2e cos (cO - ) 





cos + 8m’e' cos x; 


and therefore, integrating, 


nt + const. = 6 — 2e sin (cB - w) 
3é llm* . lime . 1s 

+ °¢ sin 2(c8 - w) - “— sin @ - —— sin f + 3me sin w, 

4 


23. We now proceed to shew that these results are not 


affected by the inclination of the moon’s orbit to the plane of 


the ecliptic. 
By Art. (14) p is the only part of R which contains ¢; 
therefore, since, by Art. (16), 


m mr 3p-1 
re ee } 


Rk - ’ 13 
r r” 2 
dR 3n'r dp 
we have age ee cE 
a pf Pa’ ) 


and, by Article (14), observing that \' = 0, 0, = @, 
p = cos d cos (0, - &), 
} do. 
‘ = 3 =—sindA “ cos (6, — 0) — cos Te sin (9, - 0’)... .(1 


a t 


Now, by Article (14), 





sin A = sin « sin (0 - v) cos « tan (6 - v), 


tan (0, - ») 








dk 
dé 


b 
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hence, differentiating and then neglecting higher powers of 
d and «, and putting 0, = 0, we find 


. —s isd 1) : 
re sin A — =csin’(@-v), cosa ee (0 — v) cos (0- v), 
& da a ' 
# therefore, by (13) and (12), 
: dp ue ; 
4 =~csin(9-,){sin (0-») cos (0-) - cos(0-») sin (0-@)}, 
Ut 
p P — sin (0 - v) sin (0' - v) cos (0 - 0), 
a 
1 3m'r — ie ; 
oe — —,— «cos (0 - 0’) sin (0 - vy) sin (0 - v)..... (14). 
d " 
Hence, if we denote by R, the part of R independent of the 
inclination, we find, integrating, 
3m'r” , oe —" 
R= R,- pyr #008 (0 - &') sin (0 - v) sin (0 - v); 
E from which it appears that the part of R depending on the 
iy inclination is of the fourth order ; therefore the only part of 
f the preceding investigation that can be affected by the in- 
i at _ ‘ d(h’ ' 
clination is the determination of h’ from — , Article (17), 
C 
and we need only look out for those terms in which the 
: , . aR . 
. : coefficient of @ is small in 76 Now 
aa ¢ 
> dR dR, 38mr , ; . ; 
t —= —1 —- —— {cos(6-0')cos(0-r)-sin(0-6')sin (0-1)} cos (0'-» 
f qd a saat lil ati 


: aR, 38m'r’ , han ; 
. | = —! —- — v.cos (20-9 - v) sin (0 - v). 
’ ; do or" 
The part multiplied by ¢ here evidently gives no term in 
which the coefficient of @ is small, we have therefore 

dR dR, 

d@ dQ’ 
and therefore the value of h*, found in Article (17), is the 
correct value when the inclination is taken into account. 


— ‘ , . dO 
Again, in finding the time, instead of h = r? —, we have, 
dt 


3 by Article (12), 
) ~ , (de <a ae 
). : h=7 \ di —2sin : dl, 
te. 4: Ge. ds 
r  -— 7) by Article (12). 
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The part involving « here is of the fourth order, and is 
therefore to be neglected altogether. 
Hence it appears that the supposition that the inclination is 


























EGER 


zero does not in any way affect the process of finding the : 
values of u and ¢ in terms of @ to the second order of small " 
quantities. It is important to remark that @ is the longitude : 


of the moon measured on the plane of her orbit ; and to in- 
troduce the longitude measured on the plane of the ecliptic, 
namely 6,, we must apply the common reduction 


6,-0=- ‘ 0@~o. 
. 4 


This explains the reason why the only term in the expression 
for the time in Airy’s Tracts, which depends on the inclina- 
tion, is the reduction just put down. 
24. We now proceed to determine the motion of the 
plane of the moon’s orbit. 
By Article (12)* we have 
dv dvr 7 dR 


a dv 
me r 4 and a = cot (O-¥) = 


therefore, by Article (23), equation (12), we have 


S = — 3m’ cos (0 — 6’) sin (0 - v) sin (0 — »), 
Si = — 3m’: sin (0 - 6’) sin (0 - ») sin (0 - v). 


In both these put 
2 sin (0 —- v) sin (6' — ») = cos (0 — 0’) — cos (0 + 0 — 2»), 


retain only the terms which rise in integration, and the 


dy 
constant term in —, and we find 


dé 


a" - an {1 - cos’ (0 - v)}, 
di 3m 


— =-— .sin 2(@ - »), 
oo eae < oo 
3m 3m 


which give y = const. - _ 6+ — sin 2(0' --»), 


3m ~— 
« = const. + = 08 2(0' - v); 





d : : 
* In Arts. (12), (8) and (9), in the expressions for = , Sin tis put by mistake 
for tan ¢ dt 
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and these values, substituted in the equation 






sin A = sin « sin (0 - »), 





















give the moon’s latitude. 


¥ 25. We shall conclude by pointing out briefly how the 

parallactic inequality and the more correct values of the mean 
motions of the perigee and node are obtained. 

In Art. (16) the term of &, involving the third power of z 


. 7 ° 3 em ‘ 
or —, 18 mI op — 3p 
r =. >. 
7" 2 
In this term neglect the eccentricities altogether, put for 
5p’ - 3p . , — , 
et its value §cos 3 (0 - 6’) + 3 cos (0 - @’), retain only 


the part multiplied by cos (@ - 6’), which rises in integration, 
and, so far as this term of R is concerned, we find 
dR 3m'r’ (8 | 
— = —— {- cos (0 - 6’) 
dr 4 {3 ( ) i he 


dR mr f_3 ' 7 ni, 
do FY 5 sin (0 OD rs 


dé 
r aR 9m a 7 
and... - = —-=-- om » = OO ~@)....5+fA) 
h? dr S aa 
dR 3 a 
or? — = - — aum’. — sin (0- @); 
do 4 “ a ( )i 
hence in ts there arises the term 
h” 
3m a 
-- —.-,cos(0- @).......... B); 
os ( ) (B) 


‘ : dr , 
and, observing that e is zero and therefore aH 7 0, in the 
‘ ( 
equation (3) there arises the term 
15m a 
(A)+(B) or - on, = eee 0’). 
8 aia 


If, therefore, F cos (0 - 0’) be the corresponding term in w, 


we have 3 
" 15 ma 
Fic -(1-my}=-— — = 
{ ( yh 8 aa’ 
og . 15m a 
» & and therefore F=e-—-—..-. 


16 a 
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Hence we have 


oy | paaaka aa 4 cos (0 - 0). 
\ 16 a J 

The term here obtained is the parallactic inequality. (See 

Airy’s Tracts.) 

26. To obtain the more correct value of the mean motion 
of the perigee, we must look out in the second member of 
equation (3) for terms of the form Ae cos (@-w). Now 
it is easy to see that, when we substitute for 7, 7’ and 0’ 
their values to the second order (see Art. 21), in order to 
expand the second member of (3) to the fourth order, 
terms of the form Ae cos (9 - w) will arise only from the 
combination of the quantities me cos (or sin) (@ - 20’ + w) and 
cos (or sin) 2(9 - @’). Hence, confining our attention only to 
the term me cos (0 ~ 20’ + w), i.e. putting for r 
15me 





f 
ayl- cos (9 — 20' + w)}, 
\ ~ 
and for 7 and 6’, a and m@ + B simply, we find (see 


Art. 17) 








d(h® 15me a ; 
Se). ~ seun’ 1-4. —* cos (0 - 26 +)! sin 2(0 - 6) 
dé 8 J 
‘ lime . 
= — 3aum’ { - a } : 
8.15. 1 
We=a See is me cos (0 - w).... $3 
an 4 ill seas 
hence, and putting au - 1 for e cos (@ - w), we find 
Fz Biekeas ——— MU oo .00. (A) 
Again, we have (see Art. 18) 
2 > 2 15 > 8 cos 2 a 
e a ak cota iF ‘}. Se ee (0-20 4 w)| 3 cos 2(0- 0) os 2(0-) 
h? dr a \ 8 J 2 
m 9.15 ; 
=— ae ~ me cos (0 - Sibkve \ 
a \ 32 J 
9.15 ; 
Sa iba ei eole — — m Pe photo el 


Lastly, we have (see Art. 19) 
dr {15me 
12 Mh 


d0 | 8 | 


REE. 


SOgENGENDE™ 


3 
; 
sa 
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1 dr dR 3m? 1l5me . eid 
— — — = — —— , —__- § - 20° +w)s -@ 
ii i ao os ee in (0 - 26 + w) sin 2 (0 - @) 


3.15 me 


and ... 


Bose eee — ——-+-— Cos (9 —w)...... 
32 a "7 
3.15, 
Biecene™ qq eee ee ee ence eee (CO), 


Hence, adding (A), (B) and (C), we have a term 
3.15 3 1 225 
pl fy +—- LI miu = —— 


4 \ 8 8 


Hence the differential equation for u (see Art. 20) becomes 


Tu ( 3m >) ae 
——-+{ 1— —]WU=....5+ ——- MU... 





m'u. 


df \ 2 16 
- du ute 3m? =. 225 ~ 
oO — - — - — errr mae 
d 2 16 
A 3m° 225 —(, 
Hence we have c?= 1 - — —- —-m’, 
2 16 
3m 225 
c=l1l- ee A 
4 32 


which determines the more accurate value of the mean motion 
of the perigee. (See Airy’s Tracts.) 
27. The more accurate value of the mean motion of the 
node is easily obtained from the equation (see Art. 24) 
dy 3m* 


SS aw em 1h eee. a 
76 : {1 — cos 2(6 - »)} 


by putting for v, in the second member, its first approximate 
value 3m 


as Tale g sin 2(0' - »,), 


which, expanding cos 2 (0' — v), and retaining only the con- 
stant part, gives 


dv 3m { ; mn .. 1 
— =—— /1- cos 2 (6 -»,) - — sin’ 2(0-»,)} 
do 4 | ar a | 
3m? 9m? 
a 32 


which determines the more accurate value of the mean motion 
of the node. (See Airy’s Tracts.) M. 0. B. 
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Il.—ON THE EQUATION (D+a)* y= X. 
Tuts short paper is intended to illustrate the formation of all 
the differential equations connected with an equation of the 
n‘* order, namely, its successive primitives of the (m— 1)", 
(n - 2)", &c. orders. One of the greatest encouragements in 
the use of the calculus of operations is the perfect view which 
it gives of all the primitives of a linear equation of the form 
(D -a)(D- b).... y = X, allowing any one of them to be 
singled out and formed independently of the rest. This is 
not so easily done when two or more of the set a, b, &c. are 
equal, the reason being that some of the primitives in such a 
case are not linear with constant coefficients. I do not intend 
to consider all the combinations of primitives which may occur 
in an equation of the general form (D+ a)"(D+by...... 
y = X, but only the particular form (D+ a)" y = X. 

To take first the simplest form, let a = 0, or y") = X. The 
discussion of this easy case at length will allow of brevity in 
treating the more general case. First, it is desirable to express 
{dz)'px in terms of {z"pxdz: let this last be P,, we have then 
the following theorems : 


I (nm + 1) fdzy" gx =2"P,-nz™'P, + &e.; 
mip p 
Ja"P,dz = P= Pons 
[ABdz = Afdx) B - A'fdzf B + A" fdzy B - &c. 


m+1- 

From the last theorem (which is an obvious extension of 
John Bernoulli’s) it follows that if A be a rational and integral 
function of the x degree, {ABdz is always integrable when 
B can be integrated m+ 1 times. Hence pz. y) is integrable 
per se, if p= A, + A+ e000 + a,.x""': its integral being 

- (n-2 mM. * 
px.y"—p'x.y” + &c. This must be called a known theorem ;* 
but it certainly is not well known. 

The mode of proceeding with y”) =X is now obvious 
enough: let our instance be y” = X, and let [Xa"dzx = P... 
Multiply successively by z’, z', z’, x’, and integrate ; we thus 
get the four primitives of the third degree, which are (the 
constants being included in the symbols P, &c.) 

m" pam > 
(1) y"=P,, 
‘ “ > 
(2) zy" -y'=P,, 
(3) 2’y" - Ixy + 2y' = P,, 
3, Ul € m y U >, ade 
(4) zy" - 382°y + 6ay' - 6y = P,. 





> 





* It is a particular case, deducible from the general method of finding 
integrating factors for linear expressions. 
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To find the six primitives of the second order, observe that 
the first side of (1) is integrable, and remains so when multi- 
plied by z or z’; also that the first side of (2), or the same 
multiplied by z, is integrable ; also that the first side of (3) is 
integrable. Here are then the processes for finding the six 
equations required, which are 

(5) y =2P,- P,, 

(6) zy'-y'=3(¢P, - P,), 

(7) xy" - 2xy' + 24y =} (2'P, - P,), 
(8) zy’ - 2y'=2zP,- P,, 

(9) 2*y' - 3zry' + 8y =4 (2P, - P,), 
(10) 2*y’ - 4ay' + Gy =2P,- P,. 

To find the four primitives of the first order, observe 
that (5) is integrable, and remains so when multiplied by z, 
and that (6) and (8) are integrable. The four equations 
required are 

(11) y =3(@’P, - 2xP, + P,), 

(12) zy’ 7 y = 5 #P, 7 4 z'P, 7 é P,, 
(18) zy’ - 2y=32°P,-42P,+}P,, 
(14) zy’ - 38y=$2°P, -2P,+} P,. 

Finally, integration of (11) gives the primitive of the order 
zero, 1 

3 2 
(15) Y= 2.3 (a i - 32 P, a 32P, ~~ P.). 


To extend this to the case of (D+ a)'y = X. Search in the 
common way for the factor which renders (D + a)'y integrable 
per se; it will be found to be «*(c,+¢7+.... +6, ,2""), 
say «“.pz. To find D". «“px (D+ a)" y, transform it into 
D" px eo Dey or D" (pxD"e“y), which, on account of the 
degree of px, appears in finite form thus, 

(p7.D"" — p'zD?+....4 p 2D) e*y 
= «* {pz(D+ay'-p'x(D+ay*+....+ pte (D+ ay} y. 
Proceed thus with (D+a)' y= X. Let fex"pardz be signified 
by P,. We have then, as before, 





n—-1 P 
T'(n+ (Dray rganen(2P, — nx" P\+n — 2 P,- -) 
, v"P —-P 
fa™ | dz on n mintl 
: ‘ m+1 


(D+ay. AB= A(D+ay'B-A'(D+ a)*B+ A(D+ayB-... 
Though it may be worth while to exhibit these theorems in 
terms of the operation D+ a, yet it is not necessary to pro- 


62 Demonstration of Dupin’s Theorem. 


ceed further ; for (D+ a)" y= X is D"(e*y) = «*X, and the 


rule of proceeding is to solve y”) = «*X, as in the manner 


ax 


above given, and then to write in the sev eral equations «"y, 
(D+a)y, «“(D+ayy, instead of y, y', y’, &c. Thus, 
ay’ — 2y' = xf Xxdzx - { Xz‘dx being one of the primitives of 
y” = X, it follows that 
xe* (D+ ay y-2e"*(D+a)y=xf Xe"xde - [Xe“z'dz, 
or ry'+(2ax- 2) y' + (@x-2a)y=e"af Xe"adz - ¢“ [ Xe"x'dz, 
is one of the primitives of 
y” + day’ + Gay + 4a°y' + a’y = X. 

The following theorem corresponds to that of integration by 
parts (D +a)" {A (D+ a) B} = AB -(D+a)' (BDA), 
which is an immediate consequence of the simple symbolical 
relation D a 

a ee 

Dra Dia 

A.D. M. 





III.—ELEMENTARY DEMONSTRATION OF DUPIN’S THEOREM. 
Ir there be three series of surfaces, such that all the surfaces 
of each series cut the surfaces of the other two series at right 
angles, the theorem to be proved is that the lines of inter- 
section of any one of the surfaces of the three series, with the 
surfaces of the two conjugate series, are its lines of curvature. 

Let O be any point in which three conjugate surfaces 
intersect, and let the rectangular axes OX, OY, OZ be 
perpendicular to the tangent planes of the three surfaces 


at O. Let P| a on, ree .@) 
Pe ey Serer rrr rer (a,) 
pac LY Sererererer rs 


be the equations of the three series ; and, when proper values 
are attached to A, A,,A,, let (a) be the surface touched by 
YOZ, (a,) by ZOX: ‘and (a,) by XOY. 


Hence, when z= 0, y=0, z= 0, we have 


df ei df ie 

dy dz 

df. df, 

—! = 0 —!=(Q@ ee ee ees e 
dz ” dz ae ) 


Y,_, 4. 
Z dy 


dz 
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Now, since the system is orthogonal, 


df, df, Ff df, df, Ab; df, _ 
dx dz dy dy "de de 
df, df df, ff . df, df 
ant ame ft. 5 a2 wf 
dx dx dy dy” dz dz 
df df, LF df, r df df, 
dx dx Yy dy dz dz 
Differentiating the first of these equations with respect to z, 
the second with respect to y, and the third with respect to z ; 
putting x, y, and z each = 0, in the result, and making use of 
equations (6), we have 


(28) (24) +(48) (2h) -» 
(£) (Zh) (2) (#)-« 


(ae) (aeae) * (aq) Cap) ~ 


the brackets denoting that in the quantities enclosed, z, y, z 
are put= 0. From these equations we conclude that 


af af, df, 7 
(=) , (4) ’ (z:) are each = 0. 


Now, since YOZ is the tangent plane to the surface (a) at 
O, we have 


z= ($) y +2 (sa) yz + (3) ‘| + &c.....(d) 


ar) 
dydz }’ 








for points in the surface adjacent to O. To determine ( 
we have, from (a), 
f 


> ..2 
dy df 
dz 


Hence, differentiating with regard to z, and putting z, y, 
and z each = 0, in the result, 


(ae). Gg) 





dydz 
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Hence (d) becomes 
z= : (=) y+ (=) 2\ + &e. 
2 \\dy’ dz} J 

Hence the planes of zy and zz contain the principal sections 
of (a) through O, and therefore the lines of intersection of 
(a,) and (a,) with (a), touch the principal sections of (a) at O. 
Now O may be any point in one of the surfaces (a), and 
therefore each of these surfaces has its lines of curvature 
traced upon it by the surfaces of the series (a,), (a,). Similarly 
it may be shewn that each surface, (a,) has its lines of curva- 
ture traced by (a,) and (@), and each surface (a,) by (a) and 
(a,), which is the theorem to be proved. 

By differentiating the first of equations (c) with regard to 
y and z respectively, and putting z, y, and z each = 0 in 
the results, and performing corresponding operations on the 
second and third, we obtain expressions for 


\ dx > \ dx’ ]” \dy}’ \ dy]? \dz]? \ de)’ 


d’ df,\ (4S.\ 
in terms of (33) , (3) , (2) , 


which lead directly to the expression for the curvatures of the 
principal sections of the three surfaces at O, given by Lamé, 
in his Memoir on Curvilinear Co-ordinates. P.Q.R. 








IV.—NOTE ON A DEFINITE MULTIPLE INTEGRAL. 


In the XVIII number of this Journal Mr. Boole pointed 
out the incorrectness of a theorem given by M. Catalan. The 
following pages contain a brief demonstration of the result to 
which he was led. Both he and M. Catalan made use of 
what is generally known as Liouville’s theorem, and thus 
perhaps rendered their analysis less simple than it would 
otherwise have been. 
Let us transform the integral 


fdz,....fdx, f(a,t,+....+42,) 


by the assumption 
G0 +..+.4,2,=au,.... (a = Ea’), 
and by (n - 1) other linear relations connecting z,.... 2, and 
u,....u,, and such that 
Sz = Se". 
Then, as is well known, dz,.... dz, is to be replaced by 
du,.... du,, and thus 


fda, ... [da faa, ++.+4,2,)= fdu,f au, fdu,... fdu,... (1). 
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Let the integrations on the first side of this equation include 
all values of the variables which do not transgress the limits 


Sz’ = A’, Se’ =B, 


where B is supposed to be greater than A. Then, as 22*=3w’, 
the corresponding limits on the second side of the equation 


= Se = A’, Se = B. 
Transform the integral in z by assuming 


z,=rcos@,, 2z,=rsin 0, cos 8,... z,=7rsin ,...sin 0 


n-l ; 
and that in «, by similar assumptions, 

u,=rcos 0, u,=rsin 0, cos O,...u,=rsin O,...sin O,,. 

I may be allowed to mention that this transformation, 
which appears to have been given for the first time by 
Mr. Boole, in the last number of the Journal, had occurred 
to me before I had seen his paper. His analysis leads at 
once to the conclusion that dz,...dz, is to be replaced by 


r sin’*O,...sin 0, drd0,...d0._,. 


This may be also proved by successive substitutions in the 
manner pointed out in the case of three variables by Mr. A. 
Smith, in the first volume of the Journal. 


Thus (1) becomes, since E2* = Sw’ = 7”, 


B 
[ rar fsin"*0,d0,...[d0, , f{r (a, cos 0, +-.-a, sin 9,...sin 8,_,)} 
“" (B 
= | r'dr {sin"*0', flar cos @',) d0’, [sin *6',d0.,... 
A fo U U (J 
J sin 0 9 n-2 fd Ri ei 

The limits for @ and @' are the same. 

That this equation may subsist for all values of A and B, 
it is necessary and sufficient that 


fsin"*0,d0, ... [d0,_, f{r(a, cos 0,+ ... a, sin 9, ... sin 8,_,)} 


= fsin"’0’, f(ar cos 0) d0’, fsin"*0',d0’, ... fd, , .... 2. (2). 


With respect to the limits of 6 and 6, it is not difficult to 
perceive that if 0,...0., are taken between the limits 0 and 7, 
@,...2, , will receive all the values of which they are capable, 
namely, all that included between - 7 and + 7; and that the 
same set of values cannot occur more than once. But in 
order that z, may vary from 0 to — 7, it is necessary to extend 
the superior limit of 0, from 7 to 27. Thus the limits of 
6, are 0 and 2m, while those of the other variables 9,...0,_, 
are Qand w. And similarly for 0’. 

‘ 
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On the second side of (2) we have the factor 


| sin" *@',d@, sin" “6'.d0', ... [ ree). 


Let 2 be odd, then we shall have 


sin” *0',.d0', [ sin"*@'.d@', 











(27 — 4)...38.1 (n- B)4.2 27 
2 on = ——. 
(nr - 3)...4.2 (n= 4)...8. n—-3 
Similarly 
eee ee ay ae 23 , 
| sin" °O'd0', | sin""@'.d0',=——,  &c. = &e. 
» 0 v0 ; 2z-ov 
Lastly dO, = 2m. 
oi u—-l n-1 
= s * Qmr * 
Thus (A) = 2 3 i m7 1 = me oe 
n—- (n — 2 — 
awe tifoos - 81.85 FY 





Again, if m be even, we get 





T 
A)=2 
( ) [nrn~ 1 ‘% (n —] : \ 31 
{——— — —— -2]}...33 
\ 2 / 2 } ‘ 


. ae v(t , . 
And this, multiplied by ro" , or unity, gives, as before, 


(A) 








and thus (2) becomes, making r = 1, 


[" sin”*6.d6, i sin"*6,d0, . of dO, f(a, cos 0, + ... &e.) 


nm—1\ 
ty 
2) 
This is the general result given at the end of Mr. Boole’s 
paper, and which includes the others. 2.LE. 


[sine *0'.d0', f(a cos 0)... .(3). 
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1V.—NOTE ON SOME POINTS IN THE THEORY OF HEAT. 


In problems relative to the motion of heat in solid bodies, the 
initial distribution, which is entirely arbitrary, is usually one 
of the data. When this is the case, and the circumstances in 
which the body is placed are known, the distribution at any 
subsequent period is fully determined, and if our analysis had 
sufficient power, would become known in every Case. The 
solution of the problem would be an expression for the tem- 
perature of any point in the body in terms of the co-ordinates 
of the point, and the times measured from the instant at 
which the distribution is given. 

It is in many cases an interesting investigation to examine 
what this expression becomes when negative values are 
assigned to the time. If, for a particular negative value 
— 7, we find that the expression gives an actual arithmetical 

value for the temperature of every point in the body, the 

distribution represented will be such that, if it had existed 
at a time r before the initial instant, the given initial dis- 
tribution would have been produced by the spontaneous 
motion of the heat 

It is clear, however, that the arbitrary initial distribution 
may be of such a nature that it cannot be the natural result 
of any previous possible distribution, or that it cannot be any 
stage except the first in a system of varying temperatures. 
This, for instance, will be the case if there be any abrupt 
transitions in the initial temperatures of adjacent points, or if 
the curves representing the initial temperatures of points 
situated along any straight line through the body, have cusps 
or angular points ; for, though we may suppose such a dis- 
tribution to be arbitrarily made, it could obviously not be 
produced by the spontaneous motion of heat from any other 
preceding distribution, and all such abrupt transitions or 
angles which may exist in an initial distribution, will dis- 
appear instantaneously after the motion has commenced. 

If, however, in any case, a complete solution of the problem 
has been obtained, that is, if an expression for the tempera- 
ture of any point, in terms of its co-ordinates, and the time 
has been found, this expression must assume some form, when 
negative values are given to the time; and it is my object 
here to examine the nature of this form, when the initial 
state is such as not to be deducible from a previous distribu- 
tion. 

The simplest case is that of the linear motion of heat, and 
as from it we are enabled to understand the nature of the 
F2 
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problem in its most general form, I shall, for the present, 
confine myself to this case. The case of a thin rod, protected 
from any lateral radiation, or of an infinite solid, in which the 
temperature is distributed in parallel isothermal planes, is 
what is usually contemplated in speaking of the linear motion 
of heat; but the case of a thin rod losing heat by radiation 
from its sides may also be readily reduced to the same as the 
two former. For simplicity, however, we may suppose that 
there is no lateral radiation. 

Let a fixed point O in the rod be considered as origin, 
let the distance to any point P in the rod be z, and let the 
temperature of P, at the time ¢ from the initial instant, be ». 
The equation of the motion of the heat will be 

dv d’v (1): 

dg ); 
if, for brevity, we choose the unit of heat such that the con- 
ducting power of the body, referred to a unit of its volume, 
may be unity. 

The complete integral of this equation is 

= DA" cos (Mx + M,). 00 eevee (2) 

Let ,v be the initial temperature of P. Then, putting 

t= 0, in (2), we have 

go = TA, cos (2 + M,) oo 000 eveee (3). 
Now, whatever be the arbitrary nature of the function ,v, 
whether continuous or discontinuous, it has been shewn by 
Fourier that it may be represented by a series such as the 
second member of (3), and he has shewn how, when the 
value of vw, corresponding to any value of z, is given, the 
constants A;, m;, ”; may be determined. The series thus 
found is an actual quantitative representation of ,v, and is 
necessarily convergent. In fact, the proof of its convergence 
in every case must be included in the demonstration of the 
possibility of repr esenting any function ,v by the series (3). 

Now the expression for v at time ¢, is “found by multiplying 
the respective terms of the second members of (3), by the 
quantities emt, gm rt ir ssiskers ibe 
a series which converges aia t is positive, (m,, m,, &c. 
being arranged in ascending order of magnitude). Hence, 
when ¢ is positive, the series for v is still more convergent 
than the series for ,v, and therefore (2) gives a convergent 
series for v for every positive value of ¢. If, however, the 
time considered be r ~. the initial instant, the series (a) 
will become ‘she ah Tree rere 
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which is divergent. Hence it will depend on the degree of 
the convergence of the series for ,v, whether the expression 
of v in this case be convergent or divergent. In the latter 
case the distribution represented will be impossible, and 
therefore there will be no distribution from which the initial 
distribution can be derived, by spontaneous motion, in the 
time r. There are three cases of the initial distribution, 
which we must specially consider. 

1. If the convergence of the expression for ,v, or of the 
series of coefficients A,, A,, &e. 


be of a lower order than that of the series (a); that is, if the 
ultimate convergence of this series be less than that of (a), for 
every positive value of ¢, then for any finite value, however 
small, of r, the series 


2 2 
A,e” *, A,<™ ", &e. 


will be ultimately divergent. Hence, in this case, for any 
time, however small, before the initial instant, the distribu- 
tion will be impossible, and therefore the given distribution 
cannot be any stage but the first in a system of varying 
temperatures. 

2. If the convergence of the series for ,v be ultimately of 
the same order as that of the series (a); that is, if the coeffi- 
cients A, , A,, &c. can be put under the forms 


2° 2 
ae mT : a,€ mT . &e., 


where a,, a,, &c. form a series which ultimately has a con- 
vergence of a lower order than that of the series (a), according 
to the definition previously given; then, at a time +’ before 
the initial instant, the distribution is represented by 


v' = La; cos (mz + n;), 


which belongs to the first class considered, and is therefore 
an essentially primitive distribution. Hence, in this case, a 
finite age, +’ may be assigned to the given initial distribution. 

3. If the convergence of the series for ,v be of a higher 
order than that of the series (a); or if, for any finite value of 
r, however great, the series 


2 2 
Ae", Age, &e. 
is ultimately convergent; then, for every finite negative value 
of ¢, the series for v will converge, and therefore v will have 


a possible expression. Hence, in this case, no limit can be 
assigned to the age of the initial distribution. 
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The simplest case of distributions which belong to the third 
class, is that in which the expression for ,v is composed of a 
finite number of terms of the form 


A cos(mz +n); 


but we may also readily form series of an infinite number of 
terms, which will be sufficiently convergent to satisfy the 
condition stated for the third class, and therefore any initial 
distribution represented by such a series may be deduced 
from distributions existing previously for any length of time, 
however great. 

These remarks are sufficient to indicate the nature of the 
question proposed. The details of the convergence or 
divergence of the series employed will depend on the values 
which must be assigned to m,, m,, &c. for satisfying the 
conditions of any particular problem. 

If, instead of the initial distribution of heat in the rod, from 
—«@ to + being given, we have only a part, that for instance 
on the positive portion of the rod, we must, to make the pro- 
blem determinate, have some other condition given. In such 
cases, the part of the rod, over which the initial distribution 
is not given may be conceived to be removed, and the second 
condition will then generally relate to the extremity of the 
part of the rod considered. In a previous paper (On the 
Linear Motion of Heat, vol. 111. p. 170), I have shewn how a 
solution of the problem may be obtained by determining the 
initial distribution which must be made on the negative part 
of the rod, in order that the condition relative to the zero 
point may be fulfilled. Thus, if the temperature at the zero 
point, or the extremity of the part of the rod originally given, 
be constrained to be a given arbitrary function of the time, 
this constraint may be effected by producing the rod inde- 
finitely in the negative direction, and impressing on the part 
produced a certain initial distribution, determined by equa- 
tions (6) and (10) of the Article already referred to, in terms 
of the arbitrary function of the time, and the given initial 
distribution on the positive side. In many cases, however, 
the distribution so determined will be impossible, either over 
the whole extent of the negative part of the rod, or from 
some point at a finite distance on the negative side, to an 
infinite distance. Notwithstanding this impossibility, the 
solution of the problem obtained by the method described 
above will still give a finite possible expression for the 
temperature of any point on the positive side, at any time 
subsequent to the initial instant, which will be the tempera- 
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ture that would be actually assumed by the point, if the 
temperature of the zero point had, by some external applica- 
tion, been constrained to satisfy the given condition ; and the 
only thing indicated by the impossibility of the distribution 
on the negative side will be that this constraint cannot be 
effected actually, by adding the negative part of the rod, with 
a certain initial distribution of heat, to the given positive 
part. It is unnecessary to enter into this question separately 
here, as the details are very similar to those given above. It 
should be remarked however that the possibility or impos- 
sibility of the distribution on the negative part will depend 
entirely on the function of the time, which expresses the 
variable temperature of the zero point, and not at all on the 
given initial distribution on the positive part; and on this 
account, if the initial distribution on the negative part be 
impossible, all its subsequent forms will generally be impos- 
sible also. 

Before leaving this subject, it may be well to notice a 
point relative to the theory of isothermal surfaces, or surfaces 
of equilibrium, which involves considerations analogous to 
those with which we have been occupied above. 

It is a known theorem that an attraction every where 
perpendicular to any given closed surface S, may be pro- 
duced by the distribution of a given quantity of matter m 
over the surface, according to a law which is in every instance 
determinate ; but in general it will be impossible to produce 
the same effect by the distribution of matter over any surface 
in the interior of S, not coinciding with it. If, for instance, 
S were the surface of a cube, or any surface containing points 
or edges, this would obviously be impossible, and it would 
probably also be impossible in the case considered by Poisson, 
in which Sis composed of two spherical surfaces. In every 
case there will be an infinite series of surfaces of equilibrium 
without S, becoming ultimately a series of spheres having the 
centre of gravity of m for their common centre, each of which 
is such that m exerts an attraction on any point in it in the 
direction of the normal. Hence, if one of these surfaces, s, 
be given, we may not only produce an attraction every 
where perpendicular to it, by matter distributed over itself, 
but by matter distributed over S, or over any surface be- 
tween S and s, enclosing the former. Hence s is analogous 
to a distribution of heat which may be produced by a pre- 
viously existing distribution. But, unless § itself have the 
particular property we are considering in s, we cannot pro- 
duce an attraction perpendicular to s by any distribution on 
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a surface within S. Thus s is analogous to a distribution of 
heat which cannot be produced by any previous distribution 
existing at a time before it greater than a certain limit r. 
Again, in some cases the series of surfaces of equilibrium, 
of which s is one, may be continued indefinitely inwards, till 
we arrive at a surface enclosing no space, as for instance 
when s is an ellipsoid, in which case the surfaces of equili- 
brium are confocal ellipsoids, and the series within s may be 
continued till we arrive at an elliptical disc. Such a case is 
analogous to that of a distribution of heat, which may be 
produced from a distribution existing an unlimited time 
before. P. Q. R. 





V.—ON THE INTENSITY OF LIGHT IN THE SHADOW OF A 
VERY SMALL CIRCULAR DISK. 


Ir has been proved experimentally by M. Arago, that the 
central point of the shadow of a very small opaque circular 
disk is sensibly as bright as if the light were not intercepted. 
Fresnel has given a general explanation of this curious result 
in the Appendix to his Memoir on the Undulatory Theory of 
Light, (Mem. de l Institut., tom. v. p. 460. Ann. 1821-22), 
and states that Poisson has proved the same thing by direct 
integration ; but he gives no hint as to where this proof can 
be found. 

The usual method of finding the illumination of a point 
appears in this instance to present some difficulty. For if 
(fig. 1) A be the centre of the disk, Q the point in the shadow 
whose illumination is to be found, and we suppose the waves 
of light plane and parallel to the plane of the disk, also if 
PQ=4u, it is easily seen ya the whole vibration at Q will 


be expressed by 
27 ofsin = 


sin — * (ut ~u).du. 


The superior limit being @, a the inferior Qa, a being a 
point on the rim of the disk. 

The integration introduces at the superior limit the rather 
doubtful expression cos #: but, supposing this to vanish, as 
the nature of the case shews that it ought, for the vibration 
at Q cannot be affected by the small waves originating at very 
distant points of the general incident plane wave, still there 
is another difficulty ; for the integration above indicated, if 


erformed, gives 
. , CX cos 2 (vt - Qa) 


for the total vibration at Q, and therefore the illumination 
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will be expressed by C*\*, which expression being inde- 
pendent of the radius of disk, is considered to prove the 
proposition. But the difficulty is, that this expression is true 
however large the disk may be, as well as when it is very 
small, which is manifestly absurd. It is true that the same 
reasons which led us to neglect the quantity cos «% would 
2a 
r 
radius of the disk is very large, and so the illumination 
would become zero, as it ought: but this does not remove 
the difficulty from the analysis ; for, according to the principle 
usually adopted, the illumination ought noé to vanish, when 
the radius of the disk is considerable ; because, since, when a 
general wave is broken up, the effect due to the small 
elementary wave is supposed to depend on the distance of 
the illuminated point from the origin of the small wave, and 
on that only, it follows that the effect of a small wave spread- 
ing from the rim of a disk whose radius is AP (fig. 2) in 
exciting a vibration at the point Q, will be the same as that 
of a wave spreading from the rim of a disk whose radius 
is AP’ in exciting a vibration at Q’, provided only that 
PQ= P'Q. In fact the common method does not take into 
account the effect of obliquity, and seems to involve the same 
error as if, in calculating the illumination of a point by a 
bright surface according to the principles of common Optics, 
we were to omit the consideration of the fact that the quantity 
of light emanating in any direction varies as the sine of the 
angle of emanation. 


teach us to neglect the expression cos — (vt — Qa), when the 


This problem has been very elegantly treated in the 
Mathematical Journal (vol. 11. p. 141), and the result there 
obtained is the same as that found above, viz. that the illumi- 
nation is measured by C*?: but the method there given 
seems to labour under precisely those difficulties which have 
been just alluded to. The fact seems to be, that the principle 
of Huyghens, which supposes the illumination of a point on 
which a luminiferous undulation is incident, to be the same 
as that which would result from the combination of spherical 
wavelets diverging with equal intensity in all directions from 
the elements of the wave-front, cannot be looked upon as a 
physical principle, but only as an artifice rendered necessary 
by the state of analysis, and which will not always represent 
the physical conditions of the problem. 


In this paper I have adopted Huyghens’ principle so far 
as this, that I consider the vibration of the ether at a given 
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point to be the sum of the vibrations which would exist if 
the wave-front were divided into elements and each of these 
elements were a centre of disturbance; but I do not suppose 
the wavelets to diverge spherically, but according to another 
law which I shall investigate. 

With this view I shall first find an integral of the partial 
differential equation, which seems applicable to the case of 
light, in the particular case im which the vibration is sym- 
metrical with respect to a certain axis, and the effect rapidly 
diminishes as we recede from that azis. 

I shall assume, as the equation applicable to the vibrations 
of the ether, 

d*p , {do d*» d*)\ 
a ee en re © 2 5 
dt’ \dv® dy’ dz J 


and let the axis of z be that about which every thing is sym- 
metrical: let z*+y’= p, then @ is a function of z, p, and ¢; also 


q'¢ = 4 q's xe 1 9 at 2 
dx’ dp* dp 
= =4 - y+ 2 ~Z 
dy” dp* * dp 
7? , {d’ 7’ ‘f 
re a pub « + 4o Sf +4 d9| wererra | 
dt rt dz dp* dp j 
Now if ¢ diminishes very rapidly as p increases ; that is, if 
d*p _. 
we may suppose p very small, the term p 72 will be neg- 
QP 
do 


ligeable as compared with 7% at all events we may 
, 


integrate upon that hypothesis, and if it should turn out that 
this is not the case, the method must be abandoned. 


, ( , 
Neglecting, then, the term p ae , our equation becomes 


dp* 
d* _{a@ d 
Os | ne 
dt’ | dz’ dp! 
_ oF 
Let @=usin — (vt - 2), 
r 
where u involves p only: differentiating and substituting, 
we have dn? . Sa du) 
-—vu=ai- utd — 


a” | x dp} ‘ 
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Shadow of a very small Circular Dish. 


we . Qa 
andg@=ae “ eal (vt - z)....(4). 


Z 
Let us now see whether a was rightly neglected: 


dp’ “ae? 
oot . 2 ot te 


“a 


. v - 
small when p is small, provided ——— 


Pap ' dp Na’ 
is not large: and this is a condition which is undoubtedly 
fulfilled ; for, though A is very small, a is very large indeed, 
and v* she cannot be very different from a’. ‘I think, there- 
fore, we are justified in considering that the expression (4) 
very nearly represents the true phy sical conditions of the 
problem: in order that it may do so, however, v must be 
greater than a, and for this I can give no phy sical reason ; 
I shall, however, assume that it is so, and proceed to apply 
for mula ( (4) to the question in hand. 

If we take an annulus of the plane incident wave, whose 
radius is Vp, the area of the annulus will be wdp, and if 
AQ (fig. 1) = f, the vibration at Q due to the waves spreading 
from the different points of this annulus, will be expressed by 


mw? yp? — a? 


ys = p . 2 
ma} e * @ sin — (vt - f) dp 
F r 
2 9 m2 y?—q? 

a ——? . I . 
—.eNr @ ” sin . (wt -f); 








where 6 = the radius of the disk. And therefore the illumi- 
nation will be expressed by 


4 Sy “pe 
2 A" 5 s._.dl e a ; 
wr \v’-a@ 


This expression diminishes very rapidly as 4 increases, 
but if 5 be small it approaches very nearly to the value 
2 





2 P = @° 


not intercepted. And this seems to agree sufficiently ac- 
curately with the experimental result. 

I will now apply the same method to express the illumina- 
tion at any point within the geometrical shadow. 

If we take the projection of the given point on the disk 
as the origin, and call its distance ‘from the centre g, and 


A 2 
a’ x (F<) » which is its value if 4 = 0, or if the light was 
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measure § from the line joining this point and the centre, 
the equation to the ee of the disk will be 


r’ sin’ 0 + (r cos 9 - gy = B’, 
r — 2rg cos A= - ¥f 
r=g cos 0 + (0 - g’ sin’ 5) 
p = 9° cos 20 + b’ + 29 cos 0 v(b° — g’ sin’ 8). 


2 2 2 
: Tua : 
If instead of x —— we put 2, for shortness, and call 
a 


the illumination at the point in question J, it will be easily 
seen that ae 
VI = 41 [e~aOdp, 


where the limits of p are the above value and w, and the 
limits of @ are 0 and z, if we double the integral, since only 
half the disk is included between 0 and 7; 


™ 
. VI - “f e* [92 cos 26+02+29 cos 6 d@ v(b2-g? sin® 9)] 
nd ¢ 


0 


= a ene [ven cos 26+29 cos 6 dé v(b®-g2 sin® 9)) 
n /0 
The integration here indicated cannot be performed gene- 
rally, but there are some particular cases which are worth 
noticing. 
If g be so small that g’ and higher powers may be neg- 
lected, the integral becomes 


| ewvon? dp = [ (1 — 2nbg cos 0) dO = x, 
0 #0 
or the illumination is the same as at the centre. 

If we take in the square of g, the integral becomes 


“1 eT 
| gras gene -| (1 - 2nbg cos 0 - ng’ cos 20 
0 0 2 
+ 2n*b*g’ cos’ 0) dé 
= 7 (1 + n°b*q’), 
2,2 
Ta 2 9 
T = —— en? (1 4 2n°b*q’). 
Hence, the centre is not so bright as an annulus taken very 
near to it. 
There is one other case which I shall consider, viz. for the 
geometrical shadow, or when g = 0; in this case the integral 


becomes om Rw 
| ern (142 cos 0) JQ = ¢ nb” [ e-2nh* cos 260 dd, 
0 
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which is less than 7, because e™ is <1, and 


fe~ cos 20 JQ) -| (1 — 2nd’ cos 20) d0=7, nearly, 
0 0 


since 6 is very small. Hence the line of the geometrical 
shadow will be less bright than the centre. Between these 
there must be at least one ring of maximum brightness, there 
may be more. 

How far all this agrees with experiment I cannot say, but 
at least this method of considering the problem appears to 
embrace the two great facts of the brightness of the centre, 


and the existence of rings in the shadow. HG 


VI.—ON THE MOTION OF THE CENTRE OF GRAVITY OF 
BROKEN BODIES. 
By W. Wauton, M.A. Trinity College. 

Tue fundamental principle adopted by Huyghens in the 
determination of the Centre of Oscillation or Agitation of 
Complex Pendulums, which forms the subject of the fourth 
part of the Horologium Oscillatorium, is enunciated in the 
following words: “Si pendulum e pluribus ponderibus com- 
positum, atque e quiete dimissum, partem quamcunque oscil- 
lationis integre confecerit, atque inde porro intelligantur 
pondera ejus singula, relicto communi vinculo, celeritates 
acquisitas sursum convertere, ac quousque possunt ascendere ; 
hoc facto, centrum gravitatis ex omnibus composite, ad 
eandem altitudinem reversum erit, quam ante inceptam os- 
cillationem obtinebat.” 

From the nature of the argument of this section of the 
Horologium Oscillatorium, it seems probable, although no 
where, I believe, expressly stated, and although in opposition 
to the direct meaning of the words which we have quoted, that 
Huyghens really supposed each of the disconnected molecules, 
after arriving at its point of zero velocity, to be there perma- 
nently retained ; so that the ‘centrum gravitatis ex omnibus 
composite ” should be considered as permanently assuming its 
original altitude from the date of the arrival of the molecule 
of greatest velocity at its maximum altitude. Such a modi- 
fication of the proposition is at any rate necessary to its truth. 
The corrected form of the principle which we have mentioned 
involves indirectly, and in a particular case, the principle of 
Vis Viva, a generalization which in fact resulted from the 
principle of Huyghens in the hands of John and Daniel 
Bernoulli. A demonstration of the principle of Huyghens, 
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in its unmodified form, has been attempted by Professor 
Whewell, in a work entitled The Mechanics of Engineering, 
p. 142, by the application of the general principle of Vis 
Viva. The reasoning however of the Professor is not 
satisfactory, in consequence of his having omitted to take 
account of the idea of time, an essential element of the 
problem. 

The object of this paper is to offer a few remarks on the 
motion of the centre of gravity of the disconnected mole- 
cules, supposing them not to be retained in their respective 
positions of zero velocity, but to be capable of descending 
spontaneously from their positions of instantaneous rest. It 
is not necessary to suppose the body to have been oscillating, 
before the dissolution of the connections of its molecules, 
about a fixed axis in the manner of a pendulum, but we may 
imagine it to have descended from rest in any conceivable 
way; the remarks we are about to make being generally 
applicable. 

Let m, m', m".... represent the masses of the molecules, 
h, h', h’.... their initial distances below the horizontal line 
through the initial position (G) of their centre of gravity; 
2, x, z"....their distances below the same line at the mo- 
ment of the dissolution of their mutual connections ; then, if 
u,u',u'.... denote the velocities of the molecules on the 
commencement of their independent motions, we have, by 
the principle of the Conservation of Vis Viva, 
mu + mu? + mu? +o... 

= 29g {m (a —h) + m'(a' —h’') + m"(a" -h') +....}35 
but, by the property of the centre of g che of any number 


of bodies, mh + mh! +m"h" +....= 9, 


and mzx+m'z' + m2" +....=(m+m'+m'+....)2; 
where Z denotes the distance of the centre of gravity (HZ) of 
the molecules at the moment of their mutual disconnection 
below the horizontal line through (G); hence 

mu? + mu? + mu +....= 2g (m+ m' +m" +....)z....(1). 
Again, let v, v', v’.... denote contemporaneous velocities of 
the molecules in their independent movements, and y, y', y" 
.....-corresponding contemporaneous altitudes above the 
horizontal line through (#7); also let A, kh’, k’.... denote 
the values of y, y',y".... at the instant of disconnection. 
We will commence with supposing the molecules to be 
constrained to move from the moment of detachment along 
smooth invariable lines, and not to be subject to mutual col- 
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lision ; then, by the theory of the motion of heavy particles, 
whether the paths of the molecules previous to disconnection 
be joined continuously with those which they afterwards 
pursue, or their motions be diverted at the instant of their 
mutual detachment without loss of velocity into different 
directions, 
mv* = mu* — 2mg (y-k), m'v® = m'u® — 2m'g (y' - k’),| (2) 
mo? =m’ — Imig (y’ — Kh’), ccc scecceeecces mo 
we have therefore 
mv’ + m'v® + m'v? +....= mu + mu? + mus... 
— 29 {m(y — hk) + m'(y' - kh’) + m' (y’ - kh) +... ef 
but, by the property of the centre of gravity of bodies, 
mk+mk' +mk +....=90, 
and =omy+m'y'+ my +....=(m+m +m +....)Y, 
where y denotes the altitude of the centre of gravity of the 
disconnected molecules at any time above the horizontal line 
through (#7); hence 
mv’ + m'v"? + m'v? +. 
= mu + mu” + mu? +... - 2g (m+ m' +m +...) Y + (8). 
From this last equation it is evident that y can never have a 
value greater than that given by the equation 
2g (m+ m'+m'+....) y= mu + mu + mu? +.... 
= 29(m+m'+m +....)z% from (1); 
and therefore it is impossible for y ever to exceed Z in mag- 
nitude. 

It remains to consider whether, or under what circum- 
stances, ¥ will ever become equal to z. Suppose y = ; 
then, from (3), we see that 

mv* + m'v? + m'v® +.... = 90, 
and therefore, contemporaneously, 
o=0, ve =0, wv =0.... 
or, from the relations (2), 

2g (y-k)=w, Wy -k)=u", AwWly'-k)=u"....; 
hence we see that, if ever y become equal to Z, every mole- 
cule must have at that instant a zero velocity, or be at the 
highest point of one extremity of its arc of oscillation. Sup- 
pose now that ¢, r, denote respectively, the time which m 
takes to move through the first and second portion of its arc of 
oscillation, the two portions being supposed to be separated 
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by the position of m at the instant of its detachment. Let 
(t', x’), (t’, x’), (t", x’) represent like times in relation to the 
molecules m’, m’, m'’.... Then, since 2, v', v.... must be 
contemporaneously equal to zero in order that y may be 
equal to Z, the condition for such equality will evidently be 
expressed by a class of relations of the form 
t+n(t+ r)=ten' (C+ r)attn (f+ 7r)=....(4), 

where n,n’, n’.... are all positive integers. The possibility 
of satisfying these relations by any finite values of n,n’, 7’... . 
as will be evident on a little reflection, will essentially depend 
upon the forms of the paths described by the molecules sub- 
sequently to their disconnection. In case the molecules were 
to commence their motions at the moment of detachment 
along cycloidal ares, starting at the lowest point of each, we 
should have ee ees ee oe ee 


and the relations (4) would become 
Qn+1=2n'+1=2n +15.... 
and therefore n=nN =N =....3 
or the centre of gravity of the disconnected molecules would 
rise at each semi-oscillation of each separate molecule into 
a position on a level with (G). For the sake of another 
illustration, suppose that the molecules are perfectly elastic 
and keep rebounding above a horizontal plane in vertical 
lines ; in this case we must put r, 7, 7... .all equal to zero, 
and the relations will become 
(n+ 1)t=(n"'+1)t=(m +1) s.... 
or, since u=gt, w=gt, u'=gt'.... 
(n +1) u=(n'4+1)u' =(m'+1)u' =... 3 
if, then, w, «', wu’... . be quantities relatively commensurable, 
it is plain that these relations may be satisfied by admissible 
values of n, »',n’....3; on the other hand, if the requisite 
commensurability do not subsist among the projectile veloci- 
ties, it is equally evident that the relations cannot be satisfied 
by any finite values of 2, 2’, n’....; in the one case, then, 
the centre of gravity of the disconnected molecules will from 
time to time, at finite intervals, arrive at its original altitude, 
while in the other case it will never attain to it in any finite 
time. It is needless to multiply examples of the impossibility 
of satisfying the formula (4) by any finite values of 2, n', 2... 
it being obviously rather a singular coincidence that the 
requisite commensurability should subsist without previous 
arrangement than a state of things of probable occurrence. 
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We have been supposing in the preceding observations 
that the disconnected molecules oscillate permanently in in- 
variable curves, and that they are free from a liability to 
mutual collision. We will now suppose the molecules, con- 
sidered as perfectly elastic, to be capable of assuming any 
irregular motions whatever, and of impinging against each 
other. Under this more general consideration of the possible 
movements, the formula (2) will be no longer applicable, the 
values of v, v’, v.... obviously depending not only upon 
u,u',u.... and the altitudes of ascent, but likewise upon 
the effects of collision. ‘The formula (3) however, as we 
know by the application of the Principle of the Conservation 
of Vis Viva, still continues to hold good; we may there- 
fore shew, as before, that it will be impossible for the centre 
of gravity of the molecules to rise to its original altitude, 
unless the velocities v, v', » ... . become simultaneously equal 
to zero. Let ¢, denote the interval from the instant of 
detachment to the arrival of the molecule m at its first 
position of zero velocity ; ¢,, ¢,, ¢,....representing the inter- 
vals between its successive arrivals at such positions: let 
t,’, t,, t,....denote similar intervals in respect to the mole- 
cule m', and so on for the rest. ‘Then, in order that the 
centre of gravity may attain to its original altitude, we must 
have 
£6666 4k +E 8 a8 OE lo, 


1 2 3 n’ 1 2 
the suffixes n, n’, x ....denoting certain integral numbers. 

The probability that for any finite values of n,m’, n’.... 
these indefinitely numerous relations should be verified, is 
evidently indefinitely small. We may therefore lay it down 
as a general practical principle, that if a mass of matter, either 
in its descent from rest or at the instant of its starting, be 
dismembered into its constituent molecules, which subse- 
quently pursue paths prescribed by the accidental circum- 
stances of gravity or collision, even supposing the molecules 
to be perfectly elastic and discarding all consideration of 
friction or atmospheric resistance, its centre of gravity will 
never again rise to a level with its original position. 

The truth of this practical principle, it will be evident, 
depends not upon any physical impossibility of the arrival of 
the centre of gravity at its original altitude, but upon the 
infinite mathematical improbability of the conspiracy of the 
necessary circumstances within any finite time. ‘The observa- 
tions which we have made above, in relation to the move- 
ments of the disconnected molecules of descending bodies, will 


G 
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evidently be equally applicable to the case of fluids flowing 
down from a higher into a lower level, the centre of gravity 
of which will accordingly oscillate in a lower basin without 
ever reaching its original altitude. 

If instead of supposing the mass to be dismembered into its 
constituent molecules, we imagine it to be broken merely into 
fragments of finite magnitude, since the equation (3) will still 
hold good, we shall still have relations such as those of (5) 
for the conditions of the rise of the centre of gravity to its 
original height ; it is important to remark, however, that the 
number of these relations will, according to the new hypo- 
thesis, be finite; for if three molecules of each fragment, 
which are not in a single straight line, ever assume zero 
velocities, the whole of each fragment will be in a position of 
instantaneous rest. ‘Thus, if there be p fragments, there will 
be only 3p—1 relations of condition ; if the mass be supposed 
never to have fallen to pieces at all, then, according to the 
same reasoning, we see that there will be only two relations 
of condition. ‘Thus we come to this general conclusion, that 
if a mass of matter in descending from any position fall to 
pieces, the improbability of its centre of gravity ever rising 
to its original altitude will be the greater, the greater be the 
number of the fragments ; and that, when the fragments 
degenerate into molecules, the improbability will amount to 
a practical impossibility. It may perhaps not be superfluous 
to observe that, if all the particles of a descending body be 
constrained to move in parallel planes, the whole body will 
be at rest if any two molecules, not in the same horizontal 
line, be at rest, and therefore the equations of condition will 
be reduced to one. If the body oscillate about a fixed axis 
as a complex pendulum, then if any one molecule be at rest 
the whole body will be for a moment without motion; and 
therefore there will be no relations of condition, or the centre 
of gravity will undoubtedly rise at the end of each oscillation 
to its original altitude. 





VII.—ON THE INVERSE CALCULUS OF DEFINITE INTEGRALS. 
By Georce Boo.e. 
Tis paper must be considered as the continuation of a 
paper on the Zransformation of Definite Integrals, (Journal, 
No. Xvil. p. 216.) 
M. Liouville has given the following theorem as the basis 
of an inverse calculus of definite integrals, viz. 
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dz x" (x + a) = (-)" r(n)( & ) p(a)....(28), 
Jo da, 
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a demonstration of which will be found in an able article in 
this Journal, vol. 1. p. 118. It is there observed, that if the 
form of @(x) be such as to render the first side infinite, the 
second may be made to agree with it by the aid of the com- 
plementary function, so that the theorem is true for all forms 
of p(x). M. Liouville, I believe, considers his result as only 
true under a limitation, and this appears to me to be the 
correct view. ‘There are at least two distinct cases of the 
theorem, each of which admits of a convenient transformation 
in the case of » being fractional. 

Let us consider the integral {dz a""p (2 +a), assuming in 
the first instance that @(z) may be developed in ascending 
powers of x, and that the first exponent exceeds - 1. 

Put x = - az, and let a = «*, then 


Ida x""o (a + x) =(-Ja" fdz2""  {a(1 - 2)} 
= (-)e"® fdzz"" {e°(1 — 2)} 


d 
= (-)'e"? fdz2"' (1 - z)™ o(e®). 
If we take for the limits of z, 0 and 1, the corresponding 
limits of z will be 0 and - a, and the above gives 


C(n)T (5 + 1) 


r(S. n + 1] 





| ‘dex (a+ 2) = (-)re  (e"). .. .(29), 


a result unquestionably true for all positive values of ». Let 
n be an integer, then 


[ia 2"'o(a+ 2) 
° — (_ yn pr d (@ . ae oom he 1) 6 
= (-)’e""I (3 { n)( ap n ) ‘\ 26 li  (e*) 


(-)' P'(n) ee 5: (5-n i)! ep (e") 
( \¢ 


' ery". 
= (-)" [(n) ‘0 (3) a"p (a) 


a \* 
-(-yYT mae \. em: 
(-) r(w)( 5) (a) ) 


a conclusion easily verified by integration, it being observed 
that the limits of the integrations relative to @ are 0 and a. 
We might obviously regard this theorem as true for fractional 
values of x; but, without thus extending our assumptions, we 
may see from the above, that M. Liouville’s formula does 
G2 


It 











‘| 84 Inverse Calculus of Definite Integrals. 


not apply when @(z) is only susceptible of development in 
} ascending positive powers of z. 
h Let us now suppose the development of @(«) a descending 
one, and the exponents negative. 
Put z = az, a= &’, as before, P 
fda 2" (a+ x) = [dz z"" (1 + 2) o (e*). .. .(81), 


a — ‘ . 
wherein 7 * of negative interpretation. Now, by a well- 
, ; 
known theorem, (Gregory’s Examples, p. 473), 


fas I (n) U(r - n) 
dz2"'(1+2z)"= ——— teens $2), 
} 5 (32) 
provided that » be positive and less than r. Hence, com- 
paring with the second member of (31) taken between the 











same limits, a 
a ri(n)l|-—-n)} 
nO . n-l ‘ 9 nQ ) \ dé } 
g dz 2" (1+ 2)" o (e) =e om a , 
/0 r ( “i a 
dé) 


whence, observing that the limits of z are also 0 and w , 


! f{ @ 
Pin)T(-—-n} 
\_ 46 / o()...(83) 





[ dz x"" p(a+ x) = <"9 


Jo rf d\ 
\ dé) 
rin) ff-2-1\ (-2-a\\ of 
g r(m) 4 76 1). .{ 76 nm)  (*) 
\ +d aa 
=(-yr <0 }( d \. @ ee | | 4 
PE (m) eye) lage tf oe 
\-n 
= (-)"T'(n) | £ GIG) wccsveicceBedsases (34). 


/ 


Hence, it appears that M. Liouville’s formula applies when 

@ (x) involves only negative exponents in its development, 
Shomal the first with its sign changed exceeds x. The 
limits of the integrations relative to a in the second member 
are+@anda. § 6. 

Finally, let us suppose x frational. 

The object of investigation being to determine the function 
under the sign of integration, when the definite integral 
is given, let us assume that within the supposed limits, 
dx a @(a+2x) = £(a), a known function, and that we wish 
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to ascertain the form of @ (a). In the case of the meoentiang 
development we have from (29), on changing a into °°, 
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r( veki 
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W (2%) = (-)"'(n) "9 o («°), 


- (-y" mmm @ 1 


8 > 
l'(n) Ta Wb (e). «22... (86). 
ris -n+1 





\d@ 


Let ¢ be the integer next below x, being 0 when » is a 
proper fraction, 


d 
en a) ny @ d \ cd 
| rs aa *7O\ a0 ~ 1) (a> 
\d 


,(@ ) a 
I -7)T(¢-n+1) 





_ ag! 1 \d0 (<8) 
\da) J P@=n +1) 5. 1) ™ 
\d0 
2 ai 023 ; 
'- “gegen (as) Jae an TC 


a’ n+l d 
= et. a. -1 1 = i-n ( . 
G-ns1) (a) [doo vo (1 - vf P(av) 


; (\ ae - oY] 0 
9a) I(n)T(t-n +1) (33) ik (1 0 “Kae. 87). 


Similarly for the descending development of (x), and for 
negative exponents, we have, by (33), 
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(Qh) me mle CM  (c*) 2... (88 
{ ete) c(n) _( a ; 4 , 
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dé 
r{ : +n | 
as nee oD 
: "6 do) ( 29) 
I'(2) r { ad \ 
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Now let ¢ be the -— r next below x, then 


a ae ee eee, 
r(- 5+n)- (- 3 1 +n-1), \~ 90 ni) — +" é 


d d\' / da.) 
on ee 
r( 5) \ J \" a) 
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d 
Ti ly 1 
we T(n) \da da V(t-n+1) r ok 
dé 


= \i+] 1 San d an” +1 | n-j-1 .\i-n a 39 
=) C(n)T(i-n+1) (3 ) bm dv v"*" (1-0) ¥(2).. of 


. — a 
since v “y()=y (<). 


We shall somewhat simplify the above formula by changing 
2 into 7-1, and in (30) and (37), x into — 2, in order to avoid 
the negative limit. Our results will then be comprised in the 
following theorems. 
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1. If, m being a positive fraction, 


| dx x"" ¢a-2x)=y~(a), 
0 

and if ¢ being the integer next above n, 7 - 1 be less than the 
first exponent in the ascending development of (a), then 


(a) = (3) Jie o*(1 — of" p(av); 


I (n)V(¢-n) \da 


but if z be a positive integer less than the same first exponent, 


1 d \" 
p (a) = Pn) (5) (a). 


a) 
2. If, x being a positive fraction, dz x""(a+2z)=W(a), 


« VU 
the exponents in the development of (a) being all negative, 
and if ¢ be the integer next above n, then 


1 d¥ df , ’ a 
Pow ee lv Ae l-~o é-n-1 : 
o(a)=(-} IP (n) (in) (55) ia}, ee ¥(¢) 


but if » be a positive integer, 


., J ay 
(a) = (-) I (n) (=) W(a). 


The application of these formule, as it presents no diffi- 
culties, so it requires no illustration. They enable us, when- 
ever y/(a) satisfies the necessary conditions, to determine ¢ (a) 
either in the form of a definite integral, or in that of a series ; 
but for the latter object it is perhaps simpler, changing x 
into ¢, to employ at once the untransformed equations (35) 
and (38). 

I apprehend that the same principles may be extended, by 
aid of the theorem of M. Dirichlet, to the inverse calculus of 
definite multiple integrals. At present, however, I have not 
leisure to pursue this inquiry, nor am I so much concerned 
to multiply results as to establish principles. 

Lincoln, Oct. 26, 1842. 


VIII.—ON A NEW SPECIES OF EQUATIONS OF DIFFERENCES. 


THE sort of equation to which the heading refers is that in 
which the degree of the equation with respect to one variable 
is defined by the other, as in u, =u , Where uw, is a 
function of the integers z and y. A problem of some interest 
leads to such an equation, and will tend to shew the sort of 
solution which may be given. 

It is well known that the number of ways in which the 
number » may be put together is 2"': thus 4 can be con- 
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structed in 2° ways; namely, 4, 31, 13, 22, 211, 121, 112, 
1111. But in this theorem orders are counted ; thus, 112 is 
considered as different from 121, and both from 211. ‘The 
corresponding problem, namely, to find in how many ways 
the number » can be constructed, different orders not being 
considered as different ways, is of much greater difficulty, 
and I am not aware of any thing having been written 
upon it. 

In the first place, if & be not less than », the number of 
ways in which ” can be put together is the same as that in 
which x + & can be composed of numbers one of which is /: 
thus 16 can be 8+.... in as many ways as 8 can be put 
together. If then wu, represent the number of ways in 
which z can be made up of y, and numbers not exceeding y, 
the problem is solved when w,,, is found: for the number of 
ways in which z can be put together is w,,.,, ..,, whatever 
may be the value of /, from 0 inclusive. 

Again, it is obvious that the number of ways in which any 
number, as 16, can be composed of numbers no higher than 
6 and not all lower, or w,,.,, is the sum of the ways of com- 
posing 10 out of numbers not exceeding 1, 2...6 Every way 
in which 10 can be so formed, gives 16 by adding 6. Hence 


“u uw u eeee ; 
2e.y a-y.l g-y. 2 ry. y 
Meg aia ee dt Mga sa i Sa 
Hence u -u =U : 
rtl, yl ay ay yl 
or u —-U =U ’ 
z.y z-l,y-l ry.y 


an equation of the y degree with respect to 2. 
This equation gives a tolerably easy mode of constructing 
a table: for all numbers up to 10 it is as follows; uw, ,, u, ,, 
and uw, ,, being always units. 
| 2 8 4246 6 7. 8 219 





1 
2 1 | 
3 1 & 4 
4 tk @ 4 l 
Bo -ge eB 4 ] 
6 ae 2 ES ee 
7 i oS 48 2 Ff Y 
8 l i » 4 a l 
9 ] boF «8 Seog B a.) 
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The rule being, after the units are written down, from the 
y number of the 2 line by adding the (y - 1)" number of 
the (x - 1) line to the (z - y) number in the z" line. ‘Thus 
10 can be formed out of numbers one of which is 4, and no 
other exceeds 4, in 9 ways: these 9 ways are 442, 4411, 433, 
4321, 48111, 4222, 42211, 421111, 4111111. 

As soon as we can express w,,_,., In terms of and y, we 
have then an equation of the y™ order for determining w,, | 
Say that wv, =V,. y» we have then, 

a — =V 


rey r-ysy y 
which, considering y as a constant, is of the form 
w-w,=V,. 
z 
From this consideration we can arrive at the general form of 
the solution, and this is all I shall here attempt. First we 
have w,,, = 1, so that 


2 x-2,2 


u —Uu = |; 
xz 


the complete integral of which, the constants being found 


from wu, ,= 0, u, ,= 1, is 
eT I . 
vw .= — + (- 1) 
kia 2 4 4 
ah , zr-1) l 1 
Therefore ‘. .-¢.,5— pel T's 
r.3 z so 9 } 4 


the complete integral of which, determining the constants 


from #, ,=0, #, ,=0, #,=1, is 
PE... Sod Sato! A 8 ("+ y") 
2,3 : ”™ —¢ ‘m4 ~— , 
la 


where (3 and y are the imaginary cube roots of unity. This 
may be expressed as follows:—The number of ways in 
which z can be composed of numbers, none greater and not 
all less than 3, is 
we 2-1 #@-4 2+3 w-4 oe 
>—-» ——, , Tr —— 
12 12 12 12 12 
according as z divided by 6 leaves a remainder 0, 1, 2, 3, 4, 
or 5. If we were now to attempt to generalise, we might 
suppose that 
u,,, is of the form A, + P,+....+P,, 


1 2 y 


b] 


where A, represents a rational and integral function of the 


n” degree, and P, is a periodic or circulating function of 


the x order, or going through a cycle of » values. But this 
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inference would not be quite correct. The terms of each 
successive integration are of the form 

definite constant x uw7X(u*V,), (u" = 1). 
When I, is a periodic term of the form Cf", B being some 
previous root of unity (say the m™), generally speaking 
wd.u"*V_ is a simple circulating function of the m kind ; 
and if V, be of the form A", the form resulting from 
integration is also 4,3”. But if =, which happens among 
the cases arising when » is a multiple of m, then u*3.u7V, 
takes the form «*A_,,, or the function of z is raised a unit in 
dimension. Hence the most general form of the solution of 
Se ee 

wi. A, + A,P, t A,,P, Tee A,,P,, 
where a,, the dimension of the multiplier of P,, is the 
integer in » - y divided by y. 

The next step gives for w,, the denominator 864, and the 

following numerator : 


62° + 182° - 27x — 39 + (27% + 27) (- 1)’ 
+ 82 (B"' + y*" - B* - y*) + 54 {(v- 17 + -v- 1}, 
this is the 144" part of 


(0) 2° + 82’, (6) 2 + 32° - 36, 

(1) 2° + 82° - 9x + 5, (7) 2° + 82° - 9245, 
(2) z+ 82° - 20, (8) 2° + 32° + 96, 

(8) 2° + 32° - 92-27, (9) a° + 82° - 9x - 27, 
(4) 2° + 82° + 82, (10) 2° + 82° - 4, 


(5) 2° + 827- 9x2--11, or (11) 2° + 82° - 92-11. 
according as x divided by 12 gives the remainder 0, 1, 2..,. 
or 11. 

It thus appears that uw,,, may be described as the in- 
teger nearest to 2° + 12, and w,, as the integer nearest to 
(x + 82°) + 144 or (x + 32° - 9x) + 144, according as @ is 
even or odd. Probably this simple species of description 
might be continued. 





IX. 





NOTES ON MAGNETISM. NO. 1. 


By R. L. Evuis, M.A. Fellow of Trinity College. 


A GEOMETRICAL construction, by means of which the action 
of a small magnet on a distant particle of free magnetism may 
be readily determined, is mentioned in a memoir by Weber, 
on the Bifilar Magnetometer (Scientific Memoirs, 11. p. 270). 
It is due to Gauss, but I do not know where he has demon- 
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strated it. A proof of it may be acceptable to some readers 
of the Journal. 

I begin by enunciating the construction in question, which 
will be easily understood without a figure. Let AB be a 
small bar magnet, ¢ its centre, P the particle of magnetism 
on which it acts. Join cP, draw PD perpendicular to it, 
meeting cP produced in D. Let cQ=}cD. Join PQ. Then 
PQ or QP (according to the sign of the magnetism of P and 
the direction of the poles of AB) is the direction of the action 
of AB on P; and oo — is its magnitude, M being the 

cP’ cQ ” } 
measure of the magnetism of AB, m that of the magnetism 
of P. 

The dimensions of the magnet being small, and its length 
in the direction of its axis being much greater than its 
breadth or thickness, we may proceed as follows. 

Conceive the magnet to be composed of a series of intense 
particles ranged along its axis. Let s be the distance of any 
one of them from ¢, yds the measure of its magnetism. Also 
let cP =r, and call the angle cP makes with cD, 0. 

Then the distance of the particle in question from P, is 
(7° + s* — 2rs cos @), and consequently its action on P is 

muds 
r+ 8° — 2rs cos 0’ 





magnetic attraction being supposed to follow the ordinary 
law. The component of this action along cD, is 


, (7 cos 8 - 8). 


(r? + s° - 2rs cos 0)° 


muds 





This is approximately equal to 
ai ak +3 = cos 0) (r cos 0 - s) ds, 
r r 
my 2 
or to —, {r cos 0+(3 cos’ 8-1) s} ds; 
r 
and therefore the total action of AB on P, parallel to cD, is 
m cos 0 3 cos’ #- 1 
—,— |uds +m ~; usds, 
r r 


the integrals being taken along ‘the ‘Whole length of the 
magnet. Consequently /yds..,.+ 








since the aggregate magnetisth ‘of’ a magnet, each element 
being taken with its ‘proper sign, is zero. Again, /usds 








| 


| 
| 
| 
| 








ne 


ee 
Fae ta 


92 Notes on Magnetism. 


being the moment of the magnetism of AB, is the measure of 
its magnetic power, or what we have called J. 
Consequently the action parallel to cD is 


Mm : 
—— (8 cos’ @-1)......... avd 
2 
The action of the element at s, perpendicular to cD, is 


muds ; 
PO _r sin 6, 


+ s° — 2rs cos 0Y° 


(7 


or, approximately, 


s ° 
+ 38 - cos 9) sin @ds. 
> 


The total action perpendicular to cD is, therefore, 
= sin @ cos 6...... eens es Oe 

The equation to the resultant of these two forces is, since the 
line passes through P, 

y-rsin@O x2-rcos8 

8 sin 0cos@ 3 cos’ @-1° 
For y = 0, or at the point Q, 

z,=r cos 0-197 sec 0 (3 cos’ 6 - 1), 
or 2z,=}7r sec 0. 
Now cD = r sec 0, therefore 
cQ = 5 cD, 
which proves the first part of the construction. 
Next, to find the magnitude of the whole action, square 

and add (1) and (2), then 


— ; t 
_ (1 - 6 cos? 8 + 9 cos* 8 + 9 sin’ 8 cos’ 6), 
F 
Mm ‘ 
or —— (1+ 8 cos’ 6} 
is the magnitude sought. Now I 
PQ =r {sin’ 0 + (cos 8 - } sec 0)°}! t 


=r(1 -24} sec’ 0) 
r(sect @ +4 245 5 or 


_i 
and cQ=I|rsec@.,... . wa hei a 
PD. .* 
eQ 


Therefore “(1 + 3 cos’ OF. 
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Consequently, if R be the magnitude of the resultant sought, 
Mm PQ 
 eP® cP’ 
which was to be proved. 

The preceding formule enable us to determine all the 
circumstances of the mutual action of two magnets, which are 
such as to fulfil the conditions of our hypothesis. 

For instance, in the memoir IJntensitas vis Magnetice 
Terrestris, Gauss has shown that if a magnet and a needle 
be placed at right angles to one another, then the moment of 
rotation of the needle due to the action of the magnet is 
approximately twice as great when the line of the axis of the 
magnet passes through the centre of the needle, as it is when 
the line of the axis of the needle passes through the centre of 
the magnet. 

In neither case has the transverse force (2) any tendency to 
produce rotation. Its effect is destroyed by the fixed centre 
of the needle. 

Let y be the distance of any element of the needle from the 
centre, mdy its magnetism. ‘Then, in the first case, we shall 

4 , , 
have @ = A approximately, (# being the distance between 
the centre of the magnet and that of the needle). Conse- 
quently the force (1) may be expressed by the formula 
Mmdy / a a ‘. . 
o_ (2 3 5 since cos’ 9 = 1 — @ nearly: 
2 at \ vu , 
(R+ y*y ‘ , 
which, neglecting y’, becomes 
2Mmdy 
we 
me ; . 2Mmydy 
the moment of this round the centre of the needle is a ‘ 
. . v 
and the total moment is, therefore, 
2MM' 
\ ae 
In the second case, # = 0 for every element of the needle, and 
the moment sought is, therefore, 
MM. Mn 
—— , since (1) then becomes - —— 
| a , 
and the sign is immaterial. ‘The moment in this case is, there- 
fore, one half of what it was before, which was to be proved. 

This result is included in the following general investiga- 

tion, in which we shall ascertain the moment of rotation due 





where M' = [mydy. 


bf 
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to the action of one magnet or needle upon another, whatever 
be their relative positions. 

The data by which we shall suppose the position of the 
magnets to be determined, are the distance between their 
centres, the angles which the axis of each respectively makes 
with the line joining their centres, and the angle between the 
axes themselves. The last element may be readily replaced 
by the dihedral angle between two planes which intersect in 
the line joining the centres of the magnets, and one of which 
passes through the axis of one of the magnets while the 
other passes through that of the other. 

Let the axis of the magnet, whose action on the other is to 
be calculated, be taken as axis of z, the centre of the magnet 
being the origin of co-ordinates. Let z, y, z be the co-ordi- 
nates of any point in the other magnet. Then, if r’=2*+y'+2’, 


2 2 a | 2 

3 cos’ §@-1= =>, and sin @ cos 8 = che! 9D 
Again, let a,b, ¢ be the co-ordinates of the centre of the 
second magnet, a, (3, y the cos angles its axis makes with the 
co-ordinate axes, p the distance of any element dp from the 
centre, X, Y, Z the forces on dp parallel to the axes of 
co-ordinates, m the intensity of dp. ‘Then 


A we 
X = im (32° - 7°) dp, 
r 





5 


‘ I 
Y= : im rydp, 
7 





o 


, soMmn 
Z=—, 


xzdp. 


Let G, H, K be the moments of these forces about lines 
drawn through the centre of the second magnet parallel to 
the axes of co-ordinates, 





3M. 
G = = x {y(z- c)- z(y - b)} dp, 
Mi —_— 
) = a {3x2 (z -— a) - (82° - 1") (2 - c)} dp, 
7 
Mi = 
K= > {(3z° - 7°) (y — b) - 8zy (x - a)} dp. 
~ 
Now r=a+ap, y=b+PBp, z=c+ yp. 
Hence, neglecting the square, &c. of p, 
3] 
(on — a (by — c3) pdp, 


R 
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M: , 
He m- {aca - (3a* - R*) y} pdp, 
_M , 
c- {(3a* - R*) B - 8aba} pdp, 
where RP=a'+ B+’. 
Integrating for p, we find, for the total moments, 














_ 38MM 
G) = = 
3MM' MM' 
(H) = — «4 (ca - ay) + Re 
ss 3MM' MM' 
(K) “2 R’ a (ap - ba a) ad a fa ’ 


M' being for the second magnet what M is for the first. 
Let L be the resultant of these three moments, then 
DT’ =(G) + (A) + (KY. 
Consequently, since a’ + (3° + y’ = 1, we shall have 


L = ou a’ {R° - (aa + 6B + eyy} 


a { MM'\ j - 
\ if } ) {aa (aa + (3 + cy) a *\ 
/MM'* 
+{ FP (8 + y'). 





Now, let 6, @' be the angles which the axes of the two mag- 

nets make respectively with the line joining their centres, 

and let » be the angle between the axes themselves. Then 
a= Rcoos @ and aa+bB+cy = Reos O, also a =cos ¢. 


Consequently 
» (MI i sales ; , 
[’ = ( =) {9 cos*@ sin*6’+ 6 cos (cos 6'cos @ - cos 9) + sin’p}, 
\ 
or L= acl {1 + 3 cos 8 —(8 cos 8 cos & - cos p)'}4. 


Let x be the dihedral angle already mentioned, then 
cos @ = cos 9 cos @' + sin @ sin @ cos x. 
Consequently the last equation becomes 
MM' 
L = - >3 
R 


the required expression. 


{1 +8 cos’ 9 — (2 cos 8 cos 6 - sin @ sin 0’ cos x) }4s 
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ADDENDUM TO ART. II. 

The method of this paper may be applied to the integration 
of any linear equation with constant coefficients, by means of 
the following theorems, the first of which answers to integra- 
tion by parts, the second to John Bernoulli’s ‘Theorem. ‘Let 
D denote the operation of finding a differential coefficient, 
then 
(D+ a)". PQ= P(D+a)'Q- (D+ a)'{P'(D+a)'Q}, 
P’ being the result of DP: and, si as D + a into 0, 

7, PQ = Pe'Q- P’e*Q+ Pe"Q.. 
which, if P be rational and integral, gives @*. PQ in terms 
of 0'Q, &e. 
Thus, to find all the preceding equations of 
(D- 1) (D - 2fy = X, 

those of the 13" degree, for instance: operate upon 

(D - 1)" {x"(D - 1) (D - 2)'y} =(D - 1)" (aX), 

(D - 2)" {2" (D - 2)*(D —- 1fy} = (D- 2)" (2"X), 
for all integer values of m aud x, from m= 0 to m= 5, and 
from x = 0 tom =7. We shall thus have the 14 equations of 
the 13 degree; and from each of these, exactly as in the 
preceding part of this paper, others of the 12", &c. degrees 
may be obtained. 

The criterion which, being satisfied, shows that 

= {P (D+ ¢ v)"y} 
will, when operated upon by the symbol (D + a)", not show 
that symbol in the result, precisely resembles the correspond- 
ing criterion when a = 0, or the case of simple integration. 

The preceding method may be equally applied to finding 
the intermediate equations of linear partial differential equa- 
tions with constant coefficients, or of equations of differences. 

With respect to equations of differences, the process 
answering to integration by parts is as follows. Let E 
signify the operation of changing z into z+ Az, then, P 
and Q uatng functions of z, 

(£ + a)'.(PQ) =P(£E+alQ- (E+ a) {AP.(E+a)EQ}, 
and the extension of John Bernoulli’s theorem is, E+ @ 
being 0, 

6" PQ) = PO'Q- AP.0°EQ+ A’Pe*E’Q-.. 
which terminates when P is rational and integral. 

In applying this to such an equation as 

(E - a)" (E - b)’....y+ X, 
instead of using z* as a multiplier, it would be more con- 


venient to use the factorial 2*'’. A. DE M. 











